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The  Naval  Ocean  Research  and  Development  Activity  has  been  heavily 
involved  with  the  design  of  the  Defense  Mapping  Agency’s  (DMA)  World 
Vector  Shoreline  (WVS)  data  product.  The  completion  of  worldwide  WVS 
coverage  was  scheduled  for  June  1989.  A  Compact  Disk  -  Read  Only  Memory 
Prototype  is  also  planned.  1 

Many  Navy  systems  will  be  using  the  WVS.  The  coastline  simplification 
cm  thinning  methods  described  will  enable  the  efficient  use  of  the  WVS  product 
in  mapping  applications  that  require  lower  levels  of  resolution.  The  distribution 
of  lower- resolution  coastline  data  sets  by  DMA  could  save  individual  method 
development  for  each  system,  ensure  that  all  systems  display  the  same  coastline 
at  the  same  map  scale,  and  greatly  speed  the  display. 


W.  B.  Moseley  .1.  B.  Tupaz,  Captain,  USN 

Technical  Director  Commanding  Officer 


Executive  Summary 


The  World  Vector  Shoreline  (WVS)  product  is  being  produced  at  a  map 
scaie  of  1:250,000.  Three  algorithms  for  simplifying  or  thinning  the  WVS  to 
lower  resolution  are  implemented,  tested,  compared,  and  evaluated.  The 
Douglas-Peucker  algorithm  is  recommended  for  generating  lower  resolution 
data  sets  for  distribution  by  the  Defense  Mapping  Agency.  Other  methods 
have  significant  speed  advantages  for  real-time  applications.  Guidelines  for 
selecting  a  simplification  algorithm  for  a  particular  application  are  presented. 
A  resolution  threshold  of  9  seconds  is  recommended  for  map  scales  from 
1:1, 000, (XX)  to  1:10,000,000.  A  resolution  threshold  of  72  seconds  is 
recommended  for  map  scales  smaller  than  1:10,000,000. 
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Simplification  Methods  for  the  World  Vector  Shoreline 


Introduction 

The  fractal  nature  of  the  coastline1  makes  it 
desirable  to  sample  the  World  Vector  Shoreline  (WVS) 
data  base  in  a  manner  related  to  the  map  scale  at 
which  it  will  be  displayed.  A  survey  of  Navy  users2 
has  determined  that  several  resolutions  of  the  WVS 
are  required. 

This  report  shows  the  relationships  among  vector 
data-base  resolution,  display  device  resolution,  and 
map  scale,  then  presents  three  simplification 
algorithms.  These  algorithms  provide  means  of 
reducing  the  data-base  resolution  to  that  which  is 
appropriate  for  a  particular  map  scale  and  display 
device  resolution.  The  discussion  is  supplemented  with 
a  series  of  maps  illustrating  the  output  of  the 
simplification  algorithms  at  several  map  scales. 

Data  Base  Resolution,  Map  Scale, 
and  Display  Device  Resolution 

Digital  electronic  display  devices  have  a  finite 
resolution.  It  is  not  possible  to  display  intormation 
in  any  finer  detail.  A  typical  display  resolution  is 
100  pixels  per  inch,  which  is  used  throughout  this 
report.  Pixel  si/e  in  geographic  seconds  is  given  by 

_.  ,  .  3,600  seconds/degree 

F'l.xel  size  =  -  » 

R  (pixels/inch)  *  S  (inches/degree) 

where  R  is  the  device  resolution  and  S  is  the  map  scale. 

Table  1  lists  pixel  sizes  for  various  map  scales.  The 
display  resolution,  /?,  is  100  pixels/inch. 

Pixel  size  does  not  necessarily  mandate  a  desirable 
resolution  for  all  mapping  applications.  It  does  offer 
an  objective  resolution  threshold,  since  any  increase 
in  data-base  resolution  is  wasted  as  the  finer  resolution 
data  collapse  into  the  pixel  resolution  upon  display  . 


Table  1.  Map  scale  vs  pixel  size. 


Map  Scale 

Inches/degree  Ratio  at  equator 

Pixel  size 
(seconds) 

16  0 

1:270.000 

2.25 

10.0 

1  432.000 

3.6 

80 

1:540.000 

4.5 

4  0 

1  1.080.000 

9 

20 

1  2.160,000 

18 

1  0 

1  4,320.000 

36 

0  5 

1  8.640,000 

72 

0  26 

1  17.280,000 

144 

0  1 

1  43,200,000 

576 

Simplification  Algorithms 

Three  simplification  algorithms  were  chosen  for 
implementation,  testing,  and  evaluation.  One  algorithm 
is  not  directly  resolution-dependent.  The  other  two 
algorithms  are  directly  resolution-dependent. 
Each  algorithm  produces  output  consisting  of  a  subset 
of  the  input  points.  The  second  and  third  algorithms 
may  be  considered  to  be  “excerpolation”  algorithms  as 
defined  by  Sharman.' 

The  first  algorithm  is  simple,  but  useful,  and  is 
known  as  the  Mh  point  algorithm.  The  first  and  last 
points  in  each  segment  are  retained.  The  intermedi¬ 
ate  points  are  selected  at  some  sampling  interval,  N. 

The  second  algorithm  was  developed  by  the  author 
and  Robert  E.  Mullen  at  the  Naval  Oceanographic 
Office  in  1985.  For  this  algorithm  a  resolution 
threshold  is  specified  in  the  units  of  the  data.  The  first 
two  points  and  the  last  point  are  retained.  The  distance 
from  each  intermediary  point  to  an  imaginary  line 
through  the  two  previously  selected  points  is  compared 
to  the  resolution  threshold.  If  the  distance  exceeds  the 
threshold,  then  the  point  is  retained.  If  the  distance 
ts  less  than  the  threshold,  then  the  point  is  rejected. 

The  third  algorithm  was  described  by  Robinson-1 
as  follows:  “The  Douglas-Peucker  algorithm  allows  the 
cartographer  to  specify  a  threshold  that  controls 
the  amount  of  simplification.  For  a  specified  line 
segment  the  two  end  points  are  connected  by  a  straight 
line  and  the  perpendicular  distances  from  all  the 
intervening  points  to  that  line  are  calculated.  If  a 
perpendicular  distance  exceeds  the  specified  threshold, 
the  point  with  the  greatest  perpendicular  distance  is 
used  as  the  new  end  point  for  the  subdivision  of  the 
original  line.  The  perpendicular  distances  from  all 
the  intervening  points  arc  calculated  and  compared 
to  the  specified  threshold.  If  at  any  time  none  of  the 
perpendiculars  exceed  the  threshold,  all  of  the  inter¬ 
vening  points  are  eliminated.  The  routine  continues 
until  all  possible  points  have  been  eliminated.” 

The  advantages  and  disadv  antages  of  the  algorithms 
are  summaiized  in  Table  2. 

The  three  algorithms  appear  in  FORTRAN  77  coded 
form  in  the  appendix. 

Map  Description 

From  Table  ),  resolutions  of  9  seconds  and  72 
seconds  were  chosen  to  test  and  evaluate  the  simplifica¬ 
tion  algorithms  and  as  likely  candidates  for  useful 
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Table  2.  Algorithm  advantages  and  disadvantages 


Algorithm 

Advantages 

Disadvantages 

Nth  point 

Uniform  sampling  rate.  Plots  look  good. 

Very  fast  and  simple. 

No  guarantee  of  spacial  resolution. 
Directionally  dependent. 

Laodrum-Mullen 

Guarantees  spatial  resolution. 

Moderately  fast.  Simple. 

Always  keeps  the  first  two  points  Variable 
point  spacing.  Directionally  dependent 

Douglas-Peucker 

A  standard  published  algorithm.  Guarantees  spatial 
resolution  with  minimum  number  of  points. 

Very  slow.  Variable  point  spacing. 

Not  simple  to  implement. 

Table  3.  Algorithm  performance  statistics. 


Algorithm 

Resolution 

Vertices 
%  of  orig. 

File  size 
%  of  orig. 

Run  time 
minutes 

Nth  point 

N  =  5 

25 

48 

42 

Nth  point 

N  =  20 

11 

42 

42 

Landrum-Mullen 

9  seconds 

29 

55 

50 

Landrum-Mullen 

72  seconds 

13 

41 

46 

Douglas-Peucker 

9  seconds 

20 

47 

76 

Douglas-Peucker 

72  seconds 

10 

40 

70 

lower-resolution  data  sets.  The  9-second  resolution 
would  be  appropriate  for  high-quality  maps  at  scales 
from  1:1,000,000  to  1:10,000,000.  The  72-second 
resolution  would  be  appropriate  for  high-quality  maps 
at  scales  less  than  1:10,000,000.  For  the  Mh  point 
algorithm,  N  was  chosen  at  5  and  20  to  produce 
approximately  equivalent  resolutions.  A  sample  WVS 
data  set  was  processed  using  the  three  simplification 
algorithms  and  the  two  chosen  resolutions  to  produce 
six  simplified  WVS  files. 

Table  3  presents  the  results  of  the  six  simplifications 
performed.  For  each  simplification  and  resolution,  the 
percentage  of  the  original  vertices  remaining  after 
simplification,  the  percentage  of  the  original  file  size, 
and  the  run  time  in  minutes  are  listed.  An  estimated 


40  minutes  of  the  run  time  was  spent  on  input  and  out¬ 
put  rather  than  simplification  processing. 

The  reduction  in  file  size  is  non-linear  because 
all  of  the  attribute  information  associated  with  the  ( 

coastline  segments  is  retained  and  no  segments  are 
completely  eliminated. 

The  full-resolution  data  file  and  the  six  simplified 
files  are  plotted  as  follows: 

Figures  1-7  0.1  inch/degree  (approx.  1:43,200.000) 

Figures  6-14  10  inch/degree  (approx.  1:4,3 20,000)  ’ 

Figures  15-21  10.0  inch/degree  (approx  1:432,000) 

Each  of  these  figure  groups  shows  an  order  of 
magnitude  zoom  on  the  upper-left  corner  of  the 
previous  figure  group. 


Map  scale  =  0.1  inch/degree  (approx. 


( 


1:43.200,000) 


I  idlin'  /.  M  I'S  at  original  resolution. 


Observations 


1 .  At  this  small  map  scale  all  of  the  simplified  coastlines  appear  as  good  or  better  than  the  original, 
even  though  the  number  of  vertices  is  reduced  to  as  little  as  10  percent  of  the  original 


i 


2 


Figure  S.  KT.S'  at  original  resolution. 

Observations 

1.  At  this  map  scale  the  9-second,  or  N  =  5,  resolution  is  appropriate  for  high-quality  maps. 

2.  I  he  72-second,  or  A  20,  resolution  may  be  adequate  for  some  applications. 
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Map  scale  =  1.0  inch/degree  (approx.  1  4,320.000) 


Observations 

1.  \t  tin'.  Luo  map  scale  the  original  resolution  is  appropriate  lor  high-quality  maps. 
1  1  lie  '/  '.count,  or  \  5,  icsolution  may  be  appropriate  t ■  .r  some  applications. 

3  1  ms  uiap  M.ale  emplU'i/es  the  differences  between  the  simplification  algoiithms. 
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Conclusions 

•  A  number  of  good  algorithms  are  available  for 
simplifying  high-resolution  vector  data  bases  as  a 
preprocessing  step  to  their  display  on  digital  systems 
at  lower  resolution. 

•  When  data  bases  are  simplified  to  the  appropriate 
resolution,  the  number  of  \ertices  is  reduced;  reduced 
storage  requirements  and  faster  display  are  the  results. 

•  Generally,  the  better  the  simplification  algorithm, 
the  more  complex  and  time  consuming  it  will  be. 

Recommendations 

•  If  processing  speed  is  most  important,  then  the 
Nth  point  and  L  andrum-Mullen  algorithms  are 
recommended, 

•  If  maintaining  a  specified  positional  accuracy 
is  most  important,  then  the  I. andrum-Mullen  and 
Douglas-Peucker  algorithms  aie  recommended. 

•  It  minimizing  the  number  of  xertices  while 
maintaining  a  specified  positional  accuracy  is  most 
important,  then  the  Douglas-Peucker  algorithm  is 
recommended. 

•  For  the  one-time  production  of  simplified  World 
Vector  Shoreline  data  bases  for  distribution  by  DMA, 
the  Douglas-Peucker  algorithm  would  probably  be 
most  desirable. 


•  All  of  the  simplification  algorithms  should  be 
made  available  to  the  WVS  users  so  that  choices  can 
be  made  for  system  optimization  and  error  budgeting. 

•  The  issue  of  feature  elimination  (the  total  elimina¬ 
tion  of  features  that  simplify  down  to  a  threshold  size) 
should  be  resolved  on  a  case-by-case  basis  depending 
on  the  application. 
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APPENDIX  A.  FORTRAN  SIMPLIFICATION  PROGRAMS 
A35$DJC1 : [ DATABASE. WVS ]NTHPT. FOR;9 


4 -APR- 1 9  8 


1 

2 

3 

4 

5 

6 
I 

8 

9 

10 

11 

12 

13 

14 
1  5 
16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


SUBROUTINE  NTHPT( xin , yin , numin , nth ,  xout , yout , numout ) 
c  simplifies  using  Nth  point  algorithm 
c  INPUTS: 

c  XIN(numin)-  real-  array  of  input  x  coordinates 

c  YIN(numin)-  real-  array  of  input  y  coordinates 

c  NUMIN-  integer-  number  of  input  points 

c  NTH-  integer-  the  point  interval  to  be  extracted 

c  OUTPUTS: 

c  XOUT(numin)-  real-  array  of  output  x  coordinates 

c  YOUT ( numi n ) -  real-  array  of  output  y  coordinates 

c  NUMOUT-  integer-  number  of  output  coordinates 

c  Jerry  Landrum,  NORDA/DMAP  (601)  688-4613,  AUTOVON  485-4613 

real  xin(numin) ,yin(numin) , xout(numin) ,yout( numin) 
i  i  n  =  1 
numout  =  0 

do  100  i = 1 , numin , nth 
i in=i in+nth 
numout=numout+l 

c  write! *,*)  '  i , numout ', i , numout 

xout(numout)=xin( i  ) 
yout(numout)=yin(i) 

100  continue 

i f ( i i n . ne . numi n )  then 

r  the  loop  ended  short  of  last  point  in  segment  so  include  it 
numout=numout+l 
xou t ( numout ) =x in ( numi n ) 
ycut(numout)=yin( numin) 

end  if 
ret  u r n 
end 
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A35$DJC1: ( DATABASE .WVS ] DIST3 . FOR ; 3  4-APR-198 


1 

REAL  FUNCTION  DIST3 ( xl , yl , x2 , y2 , x3 , y3 ) 

2 

c 

returns  the  distance  from  the  point  x3,y3  to  a  line  through 

3 

c 

the  two  points  xl,yl  and  x2,y2 

5 

dx= ( x2-xl ) 

6 

dy*=(y2-yl) 

7 

c 

vertical  and  horizontal  slopes  are  trivial  cases 

8 

i f ( dx . eq . 0 . ) then 

9 

dist3-abs( x3-xl ) 

10 

else  i f ( dy . eq . 0 . ) then 

11 

dist3*=abs(y3-yl ) 

12 

else 

13 

c 

From  analytic  geometry 

14 

c 

slope  and  intercept(b)  of  line  connecting  (xl,yl)  and 

(  x2 , y2 ) 

15 

slopel=(y2-yl)/(x2-xl) 

16 

bl=yl-( slopel*xl ) 

17 

c 

by  definition  of  perpendicular  lines,  slope  from  point 

to  line  is 

18 

slope  2«=-(l/slopel) 

19 

c 

the  intercept  of  the  perpendicular  is 

20 

b2=y3-(slope2*x3) 

21 

c 

using  substitution  we  obtain  the  point  of  intersection 

of  the  line  an 

22 

c 

its  perpendicular 

23 

x=(b2-bl)/( slope  1-slope 2  ) 

24 

y=(  (b2*slopel)-( bl *  slope 2 ) )/( slope  1-slope  2  ) 

25 

c 

the  distance  we  seek  is  the  distance  from  x3,y3  to  the 

point  of  inter 

26 

dist3  =  sqrt( ( x-x3 ) **2+(y-y3  )  **2  ) 

27 

end  i  f 

28 

return 

29 

end 

30 


i 


11 


A35$DJCl : [ DATABASE .WVS ] LANDRUM . FOR; 4 


4-APR-19 


1 

SUBROUTINE  LANDRUM ( xin , yin , numin , res ,  xout , yout , numout ) 

2 

c 

Landrum-Mul len  algorithm  for  resolution  dependent  simplification 

of  a 

3 

c 

string  of  points. 

4 

c 

INPUTS: 

5 

c 

XlN(numin)-  real-  array  of  input  x  coordinates 

6 

c 

YIN(numin)-  real-  array  of  input  y  coordinates 

7 

c 

NUMIN-  integer-  number  of  input  points 

8 

c 

RES-  real-  the  desired  resolution  in  the  units  of  the 

input 

po  i  n 

9 

c 

OUTPUTS : 

10 

c 

XOUT(numin)-  real-  array  of  output  x  coordinates 

11 

c 

YOUT(numin)-  real-  array  of  output  y  coordinates 

12 

c 

NUMOUT-  integer-  number  of  output  coordinates 

13 

c 

Jeiry  Landrum,  NORDA/DMAP  (601)  688-4613,  AUTOVON  485-4613 

1  4 

15 

c 

Examines  the  input  arrays  xin  and  yin,  dropping  out  points 

within 

di  s 

16 

c 

res  , 

1  7 

real  xin (numin) ,yin(numin) ,xout( numin) ,yout(numin) 

18 

numout  =  0 . 

19 

c 

keep  first  two  points 

2  0 

do  10  i  =  1 , mi n ( numi n , 2 ) 

21 

xout ( i ) =xi n ( i ) 

22 

yout ( i ) =yin ( i ) 

23 

numout  =  numou t+ 1 

24 

10 

conti nue 

25 

| 

26 

do  20  i--3,numin-l 

27 

d  =  di st  3 ( xout  (  numout- 1 ) , yout ( numout- 1 ) , xout ( numout ) , yout ( numout ) , 

28 

c  xin( i ) , yin( i  )  ) 

29 

i f (d . gt . res )  then 

3  0 

numout=numout+l 

31 

xout(  numout )  «=xi  n  (  i  ) 

32 

yout(numout)=yin(i) 

33 

end  if 

3  4 

I  C, 

20 

conti nue 

J  J 

3* 

r 

keep  last  point 

i f ( numi n . gt . 2 )  then 

3  f* 

numon t  =  numeu  t  + 1 

xoutfnumout l-xin(numin) 

4  o 

y  o  u  t  (  n  r.  m  o  u  t )  =  y  i  n  (  numin) 

4  2 

end  if 

return 

e  n  d 
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A350DJC1  :  I  DATABASE. WVS ] DOUGLAS. FOR; 11 


1 1 -APR- 1 ? 


1 

2 

3 

4 

5 

6 
7 
9 
9 

10 
1  1 
1  2 
1  3 

14 

15 

16 

17 

18 
iy 
20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
3  3 

34 

35 
3  6 

3  7 

38 

39 

40 

41 

4  2 

43 

44 
4  5 
4  6 
4  7 
48 

4  9 

50 

51 

52 

53 
84 

55 

56 

5  7 
58 


SUBROUTINE  DOUGLAS (xin, yin, numin, res, xout,yout,numout) 
c  Douglas-Peuke r  algorithm  for  resolution  dependent  simplification  of  a 
c  string  of  points.  Ref  Robinson , A . H ELEMENTS  OF  CARTOGRAPHY, 

c  Wiley, 1984 

c  INPUTS: 

c  XIN(numin)-  real-  array  of  input  x  coordinates 

c  YIN(numin)-  real-  array  of  input  y  coordinates 

c  NUMIN-  integer-  number  of  input  points 

c  RES-  real(for  compatibility  with  other  simplification  routines 

c  -  the  point  interval  to  be  extracted 

c  OUTPUTS : 

c  XOUT(numin)-  real-  array  of  output  x  coordinates 

c  YOUT(numin)-  real-  array  of  output  y  coordinates 

c  NUMOUT-  integer-  number  of  output  coordinates 

c  METHOD: 

c  A  stack  of  pointers  to  sections  of  the  string  requiring  further 

:  analysis  is  maintained  using  a  stack  routine.  Selected  points  ar 

flaged  for  output  in  a  local  character  array  YESNO. 

c  Jerry  Landrum,  NORDA/DMAP  (601)  688-4613,  AUTOVON  485-4613 

parameter  ( MAXLEN=10000  ) 

real  xin(numin),yin(numin),xout( numin ) , yout ( numin ) 

character*l  yesno ( MAXLEN ) 

integer  start,end 

logical  prin 

pr in= . false . 

yesno ( 1 ) = ' Y ' 

c  Initialize  to  keep  first  and  last,  can  all  others 
do  10  i=2,numin-l 

yesno ( i ) = ' N ' 

10  continue 

yesno ( numin ) = ' Y ' 

c  Begin  with  the  entire  string  of  points 
start=l 
end=numi n 

c  The  following  loop  exits  when  the  pointer  stack  is  empty,  it  should  n 
c  reach  completion. 

do  100  k=l , numin 
dmax=0 

if(prin)  write!*,*)  '  start , end sta rt , end 
c  Find  farthest  point  from  line  containing  start  and  end 
do  20  i=start+l, end-1 

l  Compute  distance  from  point  to  line  containing  start  and  end  points 

d=dist3(xin( start)  ,yin(start) ,xin(end) , y i n ( end ) , 
&  x  i  n  (  i  ) ,  yi  n  (  i  )  ) 

if  (d.gt.dmax)  then 
dmax=d 
index=i 

end  if 

20  continue 

if  ( dmax . gt . res )  then 

c  We  have  a  significant  point.  Flag  the  point,  handle  the  segments. 

yesno ( index ) = ' Y ' 

save  pointers  to  the  second  segment 

call  push (* 300 , index ) 
call  push (* 300 , end ) 
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59 

c  and  process  the  first  segment 

60 

end«index 

61 

else 

62 

63 

c  We  do  not  have  a  significant  point,  move  to  next 

segment 

S'! 

call  popt *200 , end ) 

6  5 

call  pop ( *  2 0 0 , s ta r t ) 

6  6 

end  if 

67 

60 

100  continue 

69 

200  continue 

70 

c  The  analysis  is  complete,  move  the  flaged  points 

to  output 

71 

numout=0 

7  7 

do  210  i = 1 , numi n 

7  5 

i f ( yesno ( i ) . eq . ' Y ' )  then 

74 

numout=numout+l 

75 

xout(numout)=xin( i ) 

76 

yout(numout)=yin( i ) 

7  7 

if(prin)  write!*,*)  '  output 

: '  , i , x  ou  t ( 

7  0 

end  it 

79 

2  1  0  con  i  i  nuc¬ 

OC 

8  1 

lei  u  t  n 

8  2 

300  continue 

8  3 

p r i n t ERROR  Stack  length  exceeded' 

8  4 

stop 

8  5 

86 

end 

87 

SUBROUTINE  STACK( ) 

8  8 

c  Stack  push  and  pop. 

89 

c  Call  Stack!)  to  initialize  stack  pointer  to  0. 

9  0 

c  Call  Push ( *  1 abe 1 , i t em )  to  place  an  integer  on  the 

stack 

9  1 

<:  Call  Pop  (*  labe  1  ,  i  tern )  to  pop  an  integer  from  the 

stack 

9 

c  RETURN 1  is  executed  upon  overflow  and  underflow 

. 

9  5 

save 

.  ;  * 

pa  ramot e  r ( LENSTK= 10000  ) 

i n  t  ego i  i  st  ack ( LENSTK ) 

5 > 

logical  p  r i n 

i , 

p  r i n  = . true. 

9  k 

i pt i =0 

C,  i , 

J  u  ft 

return 

1  ''  ! 

ENTRY  PUSH!*, item) 

1  '  . 

i  p  1 1  -■  i  p  1 1  + 1 

i  i .  ; 

i f ( lpt i . gt . LENSTK )  then 

,  •  4 

if(piin)  w  t  i  t  e  (  *  ,  *  )  '  STACK  OVEF.PLCW 

' 

i  0  : 

!  0 »' 

RETUKNl 

t  i  b  “ 

i st  ack (  i pt r ) - i tem 

!  C  8 

RETURN 
i  t 

\  \  ' 

’  :  i-  v  pop  (  *  ,  1 1  r-m  ? 

t  1  . 

;  f  (  i  p :  i  .  1  t  .  1  )  t  )i  e  n 

i  f  !  p  i  l  n  ) 
RETURN  1 


J  t  i  to 


.  t  ,i  '•  K 


)  t  e  m -  i  r.  t  ('  k  (  i  p  1 1  ) 
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